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Neutron resonance spectrometry (NRS) has been used to measure the temperature inside Mo 
samples during shock loading. The temperatures obtained were significantly higher than predicted 
assuming ideal hydrodynamic loading. The effect of plastic flow and non-ideal projectile behavior 
were assessed. Plastic flow was calculated self-consistently with the shock jump conditions: this is 
necessary for a rigorous estimate of the locus of shock states accessible. Plastic flow was estimated 
to contribute a temperature rise of 53 K compared with hydrodynamic flow. Simulations were 
performed of the operation of the explosively-driven projectile system used to induce the shock in 
the Mo sample. The simulations predicted that the projectile was significantly curved on impact, and 
still accelerating. The resulting spatial variations in load, including radial components of velocity, 
were predicted to increase the apparent temperature that would be deduced from the width of the 
neutron resonance by 160 K. These corrections are sufficient to reconcile the apparent temperatures 
deduced using NRS with the accepted properties of Mo, in particular its equation of state. 
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A long-standing problem with shock wave experiments 
on condensed matter is the difficulty of measuring the 
temperature of the shocked state before it is destroyed by 
release waves. Most measurements have been made us- 
ing photon emission spectroscopy (pyrometry), but many 
substances of interest to science and engineering (e.g. 
metals) are opaque in the relevant region of the spec- 
trum: infra-red through visible for shocks up to the tera- 
pascal regime. Emission from an opaque material comes 
from material within the photon skin depth of the sur- 
face, which generally cannot be maintained at the pres- 
sure of the initial shock for long enough to allow useful 
emission spectra to be collected. A transparent window 
can be placed in contact with the sample to maintain an 
elevated pressure at the surface until released from the 
free surface of the window or the projectile. However, 
the mismatch in shock impedance must be taken into 
account, along with the effect of heat conduction. 

Neutron resonance spectroscopy (NRS) has been inves- 
tigated as a fundamentally different technique for mea- 
suring the temperature inside a dynamically-loaded spec- 
imen, irrespective of its photon opacity [l| . Nuclear res- 
onances are characterized by the energy and line width 
in the rest frame of the nucleus, i.e. with respect to 
the speed of incoming neutrons relative to the nucleus. 
The resonance is manifested as the variation of attenua- 
tion with neutron energy. In an NRS measurement of a 
shocked sample, neutrons of a range of energies interact 
with a volume of material. The resonance measurement 
is a convolution of the resonance in the rest frame of a 
nucleus with the velocity distribution of the nuclei in the 
sample volume, which depends on the sample tempera- 
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ture. A pulse of neutrons of a range of energies passes 
through the specimen, chosen to have a measurable nu- 
clear resonance. The pulse has a finite duration o(200) ns 
and the spectrum varies with time, lower energy neutrons 
arriving later. The spectrum of neutrons is measured af- 
ter passing through the specimen; the temperature can 
be inferred from the width of the resonance. The reso- 
nance is also shifted in energy by the relative speed of 
the specimen with respect to the neutrons, so the spec- 
trum can also provide a measurement of material velocity 
inside the sample - a net average speed which shifts the 
centroid of the resonance. In contrast, almost all velocity 
measurements are made optically (usually by the Doppler 
shift) at the surface of the specimen. 

To perform a single-shot NRS measurement of a 
shocked state, it was necessary to collect a statistically 
significant neutron spectrum within the time for which 
the sample was in the shocked state, which is of the or- 
der of 1 /its for well-understood shocks as generated by 
the impact of disk-shaped projectiles Q. To achieve the 
necessary neutron intensity, it was necessary to design a 
dedicated 238 U spallation target and moderator, and to 
induce spallation neutrons using a pulse of 800 MeV pro- 
tons which had been accumulated in the proton storage 
ring (PSR) of Los Alamos' LANSCE accelerator. 

In any shock loading experiment, a difficulty is al- 
ways the synchronization of measurements with the shock 
event. It is particularly challenging to synchronize an 
impact event, which has a long delay during acceleration 
and coasting of the projectile, with a diagnostic pulse 
from a particle accelerator, which is generally not de- 
signed to be triggered externally with a short latency or 
high a priori precision in delivery of the diagnostic pulse. 
For the trial NRS experiments, the projectile was accel- 
erated by detonating solid chemical explosive, using the 
'Forest Flyer' design which gave an uncertainty in impact 
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time from the trigger signal of o(100) ns. The Forest Flyer 
produced a non-isotropic distribution of high speed frag- 
ments. In order to protect the LANSCE beamline from 
damage and breach of vacuum, the shock experiment was 
tilted so that the projectile and shock state in the sample 
were inclined at 55° with respect to the neutrons. The 
inclination was important in interpreting the NRS mea- 
surements, as discussed below. The separation between 
the projectile and the target was 20 mm, not 15 mm as 
reported before 

NRS experiments have been performed on the reaction 
products of detonating chemical explosive, and on Mo as 
it is a standard reference material for high pressure work 
The temperature inferred was significantly higher 
than expected from shock calculations using the best 
available equations of state for Mo [l[ . The Mo was doped 
with 1.7 at. % of 182 W, and the projectile was accelerated 
by high explosive. NRS temperature measurements were 
made on two nominally identical shock experiments. If 
the shock states were identical, these measurements could 
be combined with a root-mean-square uncertainty (prod- 
uct of the probability distribution function for each mea- 
surement), as was done previously: 883±46K. The sep- 
arate measurements differed by around the sum of their 
neutron-counting uncertainties, which is not inconsistent 
with an identical temperature. Explosively-driven ex- 
periments often exhibit some round-to-round variability. 
Variations in surface velocity of 0(10%) were observed in 
similar experiments Q, and variations in the timing of 
the neutrons with respect to the shock may lead to vari- 
ations in NRS temperature, as discussed below. Round- 
to-round variability argues for combining the two mea- 
surements as the sum of their probability distribution 
functions: 872 ± 90 K. Non-idealities in the loading are 
more likely to cause an increase than a decrease in the ap- 
parent temperature, so the lower measured temperature 
may be more accurate. The sensitivity to temperature 
on the principal shock Hugoniot of Mo to uncertainties 
in the equation of state (EOS) is believed to be too small 
to account for the observed discrepancy Q. Here we 
explain the temperature discrepancy by taking account 
of plastic heating in the sample, and by considering de- 
tails of the shock loading system used, which induced 
perturbations in the shock state that affect the apparent 
material temperature as measured by NRS. 

In an ideal planar impact experiment, the strain ap- 
plied to the sample is uniaxial. Uniaxial strain applied to 
a solid induces shear stresses, and for Mo at the ~60 GPa 
pressures of the NRS experiments the shear stresses in- 
duce plastic flow. Compared to shock temperatures in a 
material without shear strength, plastic flow causes addi- 
tional heating. The previous NRS temperatures [l| were 
calculated from the scalar EOS, and neglected the effects 
of plastic heating. 

The shock Hugoniot of Mo was calculated with and 
without the contribution from plastic work. Plastic work 
increases the thermal contribution to the EOS, generally 
increasing the pressure for a given compression. The ef- 



fect of material strength, i.e. of elastic stress and plastic 
flow, was treated self-consistently in a numerical solution 
of the Rankine-Hugoniot equations for shock compression 
[U, Q . This is necessary for a rigorous prediction of heat- 
ing. Material strength was treated using the Steinberg- 
Guinan model [6J, [7j] . The effect of plastic work was calcu- 
lated to be around 53 K at shock pressures around 63 GPa 
Q - a significant contribution, but not enough to recon- 
cile the temperature discrepancy (Fig. [T]) . This analysis 
depends on the accuracy of the Steinberg-Guinan model 
at these pressures on the Hugoniot. Supporting evidence 
is provided by surface Doppler velocimetry measurements 
made of these experiments: the onset of release from the 
peak velocity is marked by an elastic release wave of am- 
plitude consistent with the flow stress predicted using 
the Steinberg-Guinan model [|[ . The magnitude of plas- 
tic heating may be estimated simply, by multiplying the 
elastic component of stress by the change in volume as 
the sample is compressed. At 63 GPa, the compression 
on the principal shock Hugoniot of Mo is 0.85, which 
is almost entirely plastic. Mo exhibits work-hardening, 
pressure-hardening, and thermal softening, so the inte- 
grated plastic work depends on the precise deformation 
history. However, Steinberg-Guinan flow stresses Y for 
Mo generally fall in the range 1.6-2.8 GPa, following the 
convention that the elastic contribution to the normal 
stress is |Y". Thus, without accounting for the precise de- 
formation path, plastic heating should be approximately 
60-100 K. 

Pyrometry measurements have been made of the tem- 
perature of Mo on release from shocking to similar pres- 
sures, with release into a LiF window (~25 GPa residual 
pressure) and into vacuum. These temperatures were 
also higher than predicted without accounting for plastic 
heating. Pyrometry is prone to other systematic errors, 
such as thermal emission from the shocked window or 
from gas or glue compressed in the gap between the sam- 
ple and the window, and enhanced plastic heating from 
the deformation of surface features such as machining 
marks. The total power in thermal emission varies with 
the fourth power of temperature, and pyrometry mea- 
surements are often more accurate at shorter wavelengths 
where the power varies with higher powers of tempera- 
ture. Pyrometry is therefore prone to inaccuracy from 
spatial or temporal variations in temperature, whereas 
NRS measures the average temperature. However, the 
inclusion of plastic heating in Mo brought the predicted 
surface temperatures into reasonable agreement with the 
pyrometry data @]. 

The speed of the projectile was not measured in the 
NRS or pyrometry experiments, so the shock state could 
not be inferred directly from the published Hugoniot data 
for Al and Mo. The shock pressure was inferred from the 
peak free surface velocity observed in each experiment, 
using mechanical EOS p(p, e) derived from the Hugoniot 
data 0]. For a given projectile speed or for a given free 
surface velocity, the shock pressure depends on the shear 
modulus and flow stress assumed in the projectile and 
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sample. Thus the 'experimental' pressures as well as the 
predicted temperatures were adjusted when strength was 
included. 

The pressures induced by the detonation of the chem- 
ical explosive in the Forest Flyer were much higher than 
the flow stress of the Al projectile. According to contin- 
uum dynamics simulations, the design used for the Mo 
NRS experiments suffered from hydrodynamic features 
which deform and damage the projectile Q. In partic- 
ular, the case profile is likely to produce curvature of 
the projectile. The projectile would also be accelerating 
on impact. If the projectile is still accelerating, there 
must be a gradient of pressure and compression through 
it. On impact, these gradients induce a shock wave with 
a driving pressure which increases with time, leading to 
an increasing particle velocity. Unless this variation is 
taken into account, the broader peak could be attributed 
to a higher material temperature. If the projectile is not 
fiat, and particularly if there are radial variations in its 
speed at the time of impact, there will be radial vari- 
ations in pressure, temperature, compression, and par- 
ticle speed. Any radial component of particle velocity 
gives a different relative velocity compared with the neu- 
trons, going around the azimuth. This variation in rela- 
tive speed broadens the neutron attenuation peak, which 
could again be attributed to a higher material tempera- 
ture unless taken into account. Qualifying experiments 
were performed on the Forest Flyer system as used in the 
NRS experiments, but they were based on arrival time 
measurements and did not probe the detailed shape or 
density distribution of the projectile. Proton radiographs 
were subsequently obtained of a similar Forest Flyer de- 
sign, initiated with a plane-wave lens, showing curvature 
of the projectile in close agreement with the simulations 
in the region of the projectile affecting the W-doped Mo 
Q. In the NRS experiments, the explosive charge was 
initiated by 61 detonators fired simultaneously. Detona- 
tor misfires may occur, which could lead to additional 
distortion of the projectile and a higher apparent tem- 
perature; misfires are a possible explanation for the two 
different NRS temperatures. 

The simulations were used to predict the variation of 
compression, temperature, and particle velocity in the 
doped Mo as a function of time. In principle, the time- 
dependent fields could be used to simulate the neutron 
attenuation with the time-dependent neutron spectra. A 
slightly simpler procedure was adopted here, predicting 
the neutron attenuation as a function of energy at a series 
of instants in time. 

Given the spatial fields of mass density p{r), temper- 
ature T(r), and velocity u(r) at some instant of time t, 
the spectral attenuation a was predicted as 
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weight of the dopant, E r the resonance energy, 
E' = ^m n \u n {E) - u\ 2 
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where m n and u n are the mass and velocity of the neu- 
trons, and 
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where is Boltzmann's constant. 
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where E is the neutron's kinetic energy, a the natural 
cross-section, / the dopant mass fraction, A the atomic 



where u n is the direction vector of the neutron beam. 

The spectrum a(E) calculated at any instant of time 
was interpreted as an apparent temperature by fitting 
a Gaussian, from which the apparent temperature was 
calculated using Eq. [3] The calculated attenuation spec- 
trum was reproduced very well by a Gaussian, so it would 
not be possible to distinguish spatial variations by vari- 
ations in the shape of the resonance. 

Spatial variations in temperature accounted for 24 K of 
the apparent NRS temperature; spatial variations in ve- 
locity (including radial components) accounted for 124 K. 
The contributions did not combine linearly: the com- 
bined effect was 156 K. The resulting apparent temper- 
ature was consistent with the Mo NRS measurements. 
The magnitude of the pressure, temperature, and con- 
tributions to the apparent NRS temperature varied with 
time and position within the sample, because of rever- 
berations in the projectile originating with the loading 
history applied at launch. The measured NRS tempera- 
ture was corrected for projectile curvature by subtracting 
156 K. The resulting temperature was consistent with the 
calculated shock Hugoniots, and lay closer to the predic- 
tion using plastic flow (Fig. [1} . 

It can be seen that the effect of projectile flatness on 
the apparent temperature should be o(100) K in these ex- 
periments by considering the variation in the component 
of material velocity along the neutron path inclined at 
55° to the axis for a material speed in the shocked state 
of ~ lkm/s and the Forest Flyer projectile curvature of 
~ 5 x 10~ 3 /mm. The resulting shock curvature gives 
a speed variation of about ±120 m/s along the axis of 
the neutrons, which equates to a effective temperature 
for 182 W atoms of around 105 K when integrated around 
the azimuth. This simple estimate ignores the detailed 
contributions from neutron speeds and finer spatial varia- 
tions which were included above, but is of the same order 
as the rigorous simulation. 

We conclude that the NRS measurements of shock tem- 
peratures were consistent with the published EOS and 
constitutive behavior of Mo, taking into account the sen- 
sitivity of NRS to radial flow induced by the explosively- 
driven projectiles used. Plastic flow was calculated con- 
sistently with the shock jump relations. Using published 
plasticity data for Mo, plastic flow was predicted to raise 
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FIG. 1: Corrected temperature inferred from neutron reso- 
nance data, compared with the shock Hugoniot for Mo with 
and without strength. The smaller error bar is the prod- 
uct distribution (identical state preparation); the larger is the 
sum distribution (variation in states). 

the material temperature by 53 K compared with the hy- 
dro-dynamic shock Hugoniot, in the pressure regime of 
the experiments. This goes a long way to reconciling the 
measurements with temperatures expected from the var- 
ious equations of state for Mo, but does not explain the 
whole discrepancy. The Forest Flyer system as described 
is likely to exhibit significant spatial and temporal vari- 
ations in loading as applied to the sample. The overall 
effect was estimated to be around 160 K, dominated by 
the contribution from spatial variations in velocity. The 
temperatures inferred from NRS seem entirely plausible 
given the combined contribution of hydrodynamic shock, 
plastic work, and spatial variations in loading. The dif- 
ference between the apparent temperatures was similar to 
the apparent heating caused by spatial variations in the 
projectile. The difference between the two points could 
well reflect the finite reproducibility of the explosively- 
launched projectile design used. The plastic work and 
non-ideal projectile contributions to apparent tempera- 
ture were all predicted using a priori modeling, with no 
adjustments made to improve the match to experiment. 



Although the NRS technique for measuring shock tem- 
peratures has been known for some years, little further 
development has occurred because of the temperature 
discrepancy in Mo. Despite this lack of effort in devel- 
opment, the per-shot uncertainty in temperature in the 
Mo experiments was similar to that attained by mature 
pyrometry measurements. There is much scope for fu- 
ture developments of the loading system and the NRS 
diagnostic. Flatter shocks could be induced with im- 
proved explosive launchers, propellant guns, electromag- 
netic guns, or lasers. Electromagnetic and laser loading 
would allow the sample to be much closer to the neu- 
tron source, and the elimination of explosive products 
means smaller neutron losses from collision with H atoms. 
Other loading histories could be readily explored, includ- 
ing ramp and multiple-shock compression, and release 
from a shock. More recent NRS temperature measure- 
ments have incorporated a filter with resonances brack- 
eting those of the sample, e.g. Ag for 182 W, improving 
the measurement of the background during the shock ex- 
periment. NRS measurements could be made without 
doping with an element of different atomic number, pos- 
sibly using the natural isotopic composition of a material, 
subtracting the known resonance of unshocked material 
if necessary. The sensitivity of the neutron detectors can 
be improved by optimizing the thickness and composition 
of the scintillator. Temperature uncertainties of 20-30 K 
appear readily possible. It may be possible to construct 
a neutron detector giving spatial resolution. Eventually, 
the resonance signal could be precise enough to allow 
moments of the density of phonon states to be measured 
in shocked material, as has been demonstrated statically 
[l(| • Now that the previous NRS results are understood, 
the technique can be used with more confidence for the 
wide range of temperature measurements of interest in 
material dynamics. 
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